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Abstract: We describe general features of frequency-dependent charge transport near 
J> strongly interacting quantum critical points in 2+1 dimensions. The simplest descrip- 

tion using the AdS/CFT correspondence leads to a self-dual Einstein-Maxwell theory 
on AdSzt, which fixes the conductivity at a frequency-independent self-dual value. We 
describe the general structure of higher-derivative corrections to the Einstein-Maxwell 

_ theory, and compute their implications for the frequency dependence of the quantum- 

critical conductivity. We show that physical consistency conditions on the higher- 
derivative terms allow only a limited frequency dependence in the conductivity. The 

•^j frequency dependence is amenable to a physical interpretation using transport of either 

j_j particle-like or vortex-like excitations. 
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1. Introduction 

The AdS/CFT correspondence has become a powerful framework for the study of 
strongly coupled gauge theories [1, 2, 3]. While it is still in a nascent stage, an 
'AdS /Condensed Matter' duality is also being developed. That is, the AdS/CFT cor- 
respondence is proving to be a useful tool to study a range of physical phenomena 
which bear strong similarity to those at strongly coupled critical points in condensed 
matter systems. A variety of holographic models displaying interesting properties, in- 
cluding superfluidity, superconductivity and Hall conductivity, have now been studied 
[4]. Further interesting models of various types of nonrelativistic CFT's have also been 
constructed [5]. 

One advantage of the AdS/CFT correspondence is the 'uniformity' of the holo- 
graphic approach, i.e., a single set of calculations can describe the system in different 
disparate regimes (e.g., u/T — > versus T/uj — > 0). This can be contrasted with more 
conventional field theory analysis of conformal fixed points [6]. However, a surprising 
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result of the original transport calculations [7] was that the frequency dependence was 
rather trivial. In particular, the conductivity (at zero momentum) showed no frequency 
dependence, i.e., it was a constant. The authors of [7] traced the origin of this remark- 
able result to the electromagnetic (EM) self-duality of the bulk Einstein-Maxwell theory 
in four dimensions. Again this holographic result stands in contrast with those from 
more conventional field theory analysis [6, 8]. 

One perspective on these results is regard them as predictions of the AdS/CFT 
analysis on the behavior of nearly perfect fluids. Such fluids are strongly interacting 
quantum systems, found near scale-invariant quantum critical points, which respond 
to local perturbations by relaxing back to local equilibrium in a time of order H/ (/cgT), 
which is the shortest possible [6]. They are expected to have a shear viscosity, rj, of order 
r\ ~ hs/ks [9], where s is the entropy density, and many experimental systems behave 
in this manner [10]. At the same footing, we can then predict that 2+1 dimensional 
quantum critical systems with a conserved charge should have a conductivity which 
is nearly frequency-independent. Furthermore, in paired electron systems where the 
Cooper pair charge is 2e, the self-dual value of the conductivity is [11] 4e 2 /h, and this 
is close to the value observed in numerous experimental systems [12]. There has been 
no previous rationale why self-duality should be realized in these experiments, and the 
AdS / CFT theory of perfect fluids offers a potential explanation. 

Measurements of the frequency dependence of the quantum critical conductivity 
in two spatial dimensions have so far been rather limited [13, 14]. Engel et al. [13] 
performed microwave measurements at the critical point between two quantum Hall 
plateaus. Their results at the critical point do not show appreciable u dependence as 
hu is scanned through fc^T. However, they did not pay particular attention to the 
value of the quantum critical conductivity (they focused mainly on the width of the 
conductivity peak between the plateaus), and it would be useful to revisit this more 
carefully in future measurements. In any case, if confirmed, the AdS/CFT perspective 
appears to be the natural explanation for this weak frequency dependence. Graphene 
also has characteristics of a quantum-critical system with moderately strong interac- 
tions [15], and its conductivity has been measured [16, 17] in the optical regime where 
u ^> T; a frequency-independent conductivity was found, equal to that of free Dirac 
fermions. This is as expected, because the Coulomb interactions are marginally irrele- 
vant in graphene [15]. However, for oj ~ T, the interactions are expected to be more 
important, and graphene may well behave like a nearly perfect fluid [18]. A test of this 
hypothesis would be provided by measurements of the conductivity of graphene in this 
frequency regime, under conditions in which the electron-electron scattering dominates 
over disorder-induced scattering. There have also been discussions of duality in non- 
linear transport near quantum critical points [19, 20, 21]. Again, there is no natural 
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basis for this in the microscopic theory, while it can emerge easily from an AdS/CFT 
analysis [22, 23]. 

Given these motivations, it is clearly useful to understand the robustness of the 
AdS/CFT self-duality beyond the classical Einstein-Maxwell theory on AdSzi. As was 
pointed out in [7] , in many constructions emerging from string theory, the Maxwell field 
would have an effective coupling depending on a scalar field and the EM self-duality 
would be lost if the scalar had a nontrivial profile. From the perspective of the holo- 
graphic CFT, one would be extending the theory by introducing a new scalar operator, 
and couplings between the new operator and the original currents holographically dual 
to the Maxwell field. Further, the nontrivial scalar profile would indicate that one is 
now studying physics away from the critical point as (the expectation value of) the 
scalar operator will introduce a definite scale into the problem. 

However, we wish to understand the limitations of self-duality, while remaining at 
the critical point. For this, a possible approach is to simply modify the CFT through 
introducing new higher derivative interactions in the bulk action for the metric and 
gauge field, e.g., see [24, 25]. The latter are readily seen to change the n-point functions 
of current and the stress tensor in the CFT. While conformal symmetry imposes rigid 
constraints on the two- and three-point functions of these operators, they are only 
determined up to a finite number of constant parameters, e.g., the central charges, 
which characterize the particular fixed point theory [26] . These parameters are reflected 
in the appearance of dimensionless couplings in the bulk gravitational theory. Hence, 
to explore the full parameter space of the holographic CFT's, one must go beyond 
studying the Einstein-Maxwell theory and begin to investigate the effect of higher 
derivative interactions in the bulk action. This is the approach which we examine in 
the present paper. In particular, we investigate the effects on the charge transport 
properties of the holographic CFT resulting from adding a particular bulk interaction 
coupling the gauge field to the spacetime curvature - see eq. (2.6). 

Our main results for the frequency dependence of the conductivity without self- 
duality are given in Fig. 1. Here 7 is the sole parameter controlling the pertinent higher 
derivative terms in the bulk action; we will argue that physical consistency conditions 
imply the constraint |7| < 1/12. 

For 7 > 0, the frequency dependence has the same non-monotonic form as that 
expected by extrapolation from the weak-coupling Boltzmann analysis [6]: a collision- 
dominated Drude peak at small u, which is then smoothly connected to the collisionless 
^-independent conductivity at large u. This similarity implies that a description of 
transport in terms of collisions of charged particles is a reasonable starting point for 
7 > 0. 

On the other hand, for 7 < 0, we observe that it is the inverse of the conductivity, 
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Figure 1: The (dimensionless) conductivity a = g\a is plotted versus the (dimensionless) 
frequency to = uj/(AitT) for various values of 7 (the coupling (74 is defined in Section 2). 
Various consistency conditions imply that 7 G [—1/12,1/12] - see discussion surrounding 
eq. (5.14). 

i.e., the resistivity, which has a Drude-like peak at small u. Under particle- vortex 
duality, the resistivity of the particles maps onto the conductivity of the vortices [11], 
as we will review here in Section 6.1. Thus, for 7 < 0, we conclude that a better 
description of charge transport is provided by considering the motion and collisions of 
vortices. In other words, for 7 < 0, it is the excitations of the dual holographic CFT, 
obtained under the EM duality of the bulk theory, which provide a Boltzmann-like 
interpretation of the frequency dependence of the conductivity. 

An outline of the rest paper is as follows: In section 2, we review some basic 
background material, mainly to motivate the introduction of the higher derivative in- 
teraction for the gauge fields. In section 3, we calculate the charge diffusion constant 
and susceptibility for the dual CFT. We turn to the conductivity in section 4 and in 
particular, we demonstrate that in the modified theory, the conductivity is a nontriv- 
ial function of u/T. In section 5, we derive constraints that arise on the coupling to 
the new gauge field interaction by imposing certain consistency conditions in the dual 
CFT. We examine electromagnetic duality in the modified gauge theory in section 6. 
We conclude with a brief discussion of our results and future directions in section 7. 
A discussion of the Green's functions at finite frequency and finite momentum is pre- 
sented in appendix A. In particular, we examine the relationship between the Green's 
functions in the two boundary theories related by EM duality in the bulk. 
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2. Preliminaries 



As with many of the recent excursions in the AdS/CMT, our starting point is the 
standard Einstein-Maxwell theory (with a negative cosmological constant) in four di- 
mensions. Hence the action may be written as 



I = I d 4 x 



2P P \ R+ I?) 4gl FabFab 



(2.1) 



The four- dimensional AdS vacuum solution of the above theory corresponds to the 
vacuum of the dual three-dimensional CFT. Of course, the theory also has (neutral) 
planar AdS black hole solutions: 

ds " = J-2 ("/to dt2 + dx 2 + dy 2 ) + 4^ , (2.2) 

where f(r) = 1 — r^/r 3 . In these coordinates, the asymptotic boundary is at r — > oo and 
the event horizon, at r = r . This solution is dual to the boundary CFT at temperature 
T, where the temperature is given by the Hawking temperature of the black hole 

~ 3r " (2.3) 



4ttL 2 



At a certain point in the following analysis, it will also be convenient to work with 
a new radial coordinate: u = r /r. In this coordinate system, the black hole metric 
becomes 

ds 2 = y^a (-/(«) dt 2 + dx 2 + dy*) + , (2.4) 

where f(u) = 1 — u 3 . Now the asymptotic boundary is at u = and horizon at u — 1. 

As discussed in the introduction, we wish to extend the bulk theory by adding 
higher derivative interactions. As usual in quantum field theory, it is natural to organize 
the interactions by their dimension or alternatively by the number of derivatives. The 
Einstein-Maxwell action (2.1) contains all covariant terms up to two derivatives, which 
preserve parity, i.e., which are constructed without using the totally antisymmetric e 
tensor. Hence it is natural to next consider the possible interactions at fourth order in 
derivatives [27]. In all, one can construct 15 covariant parity-conserving terms using 
the metric curvature, the gauge field strength and their derivatives [27]. However, 
using integration by parts, 1 as well as the identities V[ a F bc ] = = R[ a bc]d, the general 



1 Note that we also treat the four-dimensional Euler density, R a bcdR ahcd — 4i? a f,i? ab + R 2 , as trivial 
since it does not effect the equations of motion. 
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four-derivative action can be reduced to eight independent terms 

h = J d 4 x^ [«!i? 2 + a 2 R ab R ab + « 3 (F 2 ) 2 + a 4 F 4 + a 5 V a F ab V c F c b 

+a 6 R abcd F ab F cd + a 7 R ab F ac F b c + a 8 RF 2 } (2.5) 

where F 2 = F ab F ab , F A = F a b F b c F°dF d a and the ot\ are some unspecified coupling 
constants. 

In a string theory context, one might expect all of these interactions to emerge 
in the low-energy effective action as quantum {i.e., string-loop) or a 1 corrections to 
the two-derivative supergravity action - see, for example, [28]. In such a context, 
these terms would be part of a perturbative expansion where the contribution of the 
higher order terms is suppressed by powers of, e.g., the ratio of the string scale over 
the curvature scale. From the perspective of the dual conformal gauge theory, these 
contributions would represent corrections suppressed by inverse powers of the 't Hooft 
coupling and/or the number of colours. Within this perturbative framework, one is also 
free to use field redefinitions to simplify the general bulk action (2.5). In the present 
case, field redefinitions can be used to set to zero all of the couplings except three, e.g., 
«3, «4 and a6 [27]. Examining the remaining three terms, the 0:3 and 04 terms involve 
four powers of the field strength and so would not modify the conductivity, at least 
if we study the latter at zero density. Hence we are left to consider only the o.q term 
which couples two powers of the field strength to the spacetime curvature. The latter 
will certainly modify the charge transport properties of the CFT and, as we discuss in 
detail in section 6, it also ruins the EM self-duality of the bulk Maxwell theory. 

While these string theory considerations naturally lead us to focus our attention on 
a single new four-derivative interaction, they are limited to the perturbative framework 
described above. However, we would also like to extend our analysis to the case where 
the new interactions are making finite modifications of the transport properties. In this 
case, we should think of the holographic theory as a toy model whose behaviour might 
be indicative of that of a complete string theory model. Recently the utility of this 
approach has been shown in holographic investigations with various higher curvature 
gravity theories - see, for example, [25, 29, 30, 31, 32]. Further, while the couplings 
of the higher derivative interactions are finite in this approach, consistency of the dual 
CFT prevents these couplings of from becoming very large, at least in simple models, 
as we discuss in section 5. 

So given this perspective of constructing a toy model with finite couplings, let 
us re-examine each of the terms in the general action (2.5). The first two terms are 
curvature-squared interactions which do not involve the gauge field. Hence from the 
CFT perspective, these terms would only modify the n-point functions of the stress 
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tensor and so are not relevant to the charge transport. Again, the third and fourth terms 
involve four field strengths and so these would only modify the four-point correlator 
of the dual current. Hence, as noted above, these terms will again be irrelevant to the 
charge transport, if we limit ourselves to the case of a vanishing chemical potential. 
Considering next the term, we note that it contains two powers of the field strength 
and so will modify the charge transport. However, this term produces higher derivative 
equations of motion for the gauge field and so, as explained in detail in [33], the dual 
CFT will contain nonunitary operators. Hence we discard this term in the analysis at 
finite coupling to avoid this problem. Finally, the last two terms in the action (2.5) 
also involve F 2 and again modify the charge transport. However, as we discuss in 
more detail in section 7, they only do so in a trivial way by renormalizing the overall 
coefficient of the Maxwell term. Therefore we are again naturally led to consider the 
a>6 interaction alone in studying the transport properties of dual CFT. 

Hence we will study the holographic transport properties with the following effec- 
tive action for bulk Maxwell field: 



where we have formulated the extra four-derivative interaction in terms of the Weyl 
tensor C a bcd- That is, it is constructed as a particular linear combination of the «6,7,8 
terms in the general action (2.5). This particular interaction has the advantage that 
it leaves the charge transport at zero temperature unchanged since the Weyl curvature 
vanishes in the AdS geometry. Further the factor of L 2 was introduced above so that the 
coupling 7 is dimensionless. From this action, we find the generalized vector equations 
of motion: 



Note that the AdS vacuum and (neutral) planar black hole solution (2.2) are still 
solutions of the modified metric equations produced by the new action. 

In closing this discussion, we must note that the four-derivative interaction in 
eq. (2.6) has also appeared in previous holographic studies [24, 34, 35]. In particular, 
[24, 34] considered the restrictions that must be imposed on the coupling 7 in order 
that the dual CFT is physically consistent. While [35] focused primarily on a five- 
dimensional bulk theory, there is considerable overlap between the latter and the present 
paper. In particular, [35] considered the charge diffusion constant and (zero-frequency) 
conductivity, as in section 3, and bounds arising from requiring micro-causality of the 
dual CFT, as in section 5. 




(2.6) 



V a [F ab - ^L 2 C abcd F cd ] = . 



(2.7) 



7 



3. Diffusion Constant and Susceptibility 



In this section, we calculate the charge diffusion constant and susceptibility, two quan- 
tities which control the two point Green's function of the dual current in the limit of 
low frequency and long wavelength [7]. We follow [30, 35] to extend the analysis of [36] 
to accommodate our modified Maxwell action (2.6). We begin by writing a generalized 
action which is quadratic in the field strength: 

d A xV=~9 (-l_F ah X abcd F cd \ , (3.1) 



8gi 

where the background tensor X abcd necessarily has the following symmetries, 

j^abcd y£ ab\ _cd j^cdab 2^ 

The standard Maxwell theory would be recovered by setting 

X ab cd = I ab cd = 5 a c 5 b d - 6 a % c , (3.3) 

where we can think of / as the identity matrix acting in the space of two-forms (or 
anti-symmetric matrices). That is, given an arbitrary two-form f ab = —f ba , then 
fab = \l a b cd fcd- With the generalized action in eq. (3.1), the theory of interest (2.6) is 
constructed by setting 

X ab cd = I ab cd -^L 2 C ab cd . (3.4) 

Extending the discussion of the membrane paradigm in [36] to this generalized 
framework is straightforward [30]. One defines the stretched horizon at r = r# (with 
th > ro and r# — tq <C ro) and the natural conserved current to consider is then 

J a = \n b X abcd F cd \ r=rH , (3.5) 

where n a is an outward-pointing radial unit vector. Then following the analysis in [36] , 
one arrives at the following expression for the charge diffusion constant [30] : 2 

dr 



D = -J^dV-X xtxt X 1 



r=r 



-gX 



trtr 



(3.6) 



Further applying Ohm's law on stretched horizon, the conductivity at zero frequency 
is given by [35] 



a = a(u = 0, k = 0) = \ y/^g~V-X xtxt X xrxr \ _ . (3.7) 

9a 



I r=r 



2 As noted in [30], there are two conditions required for the following general formulae to hold. The 
tensor X a i, cd is: i) nonsingular on the horizon and ii) 'diagonal' in the sense discussed in section 6. Of 
course, in the present case, both of these requirements are satisfied by eq. (3.4). 
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Lastly, the susceptibility is easily determined using the Einstein relation D = <Jq/x- 
Combining this relation with eqs. (3.6) and (3.7), an expression for x is easily read off 
as [35] 

dr 



X 



-i 



-9l 



-gX 



trtr 



(3i 



Of course, if one replaces X a b cd = I a b cd as in eq. (3.3), then these expressions reduce to 
the expected results for Einstein-Maxwell theory, e.g., see [7]. 

In the present case, we are interested in X as given in eq. (3.4) where the Weyl 
tensor is evaluated for the planar AdS black hole (2.2). Hence we find 

] y /-X xtxt X xrxr \_ =1 + 47 and 



lr=r 



L 2 r 



-gX 



trtr 



^o7 



Combining these expressions in eq. (3.6), we find the diffusion constant to be 



D 



47 



I671T7V3 



arctan 



7 



1/3 



+ log 



I-87 



(1 _ 2 7 V3)3 



(3.9) 



(3.10) 



v/3 7 V3 

A plot of this result is given in Fig. 2. If we consider 7< 1, this expression simplifies 
to 

Z)-^(l + 67 + ^V + 0(7 3 )) . (3.11) 

A perturbative result for D to linear order in 7 was presented in [35] for arbitrary 
dimensions and our results above match that for the case of a three-dimensional CFT. 



Next using (3.7), we find 



0o 



9l 



47). 



(3.12) 



Note that this expression is the exact result for arbitrary 7. The simple 7-dependence 
appearing in the conductivity contrasts with the complicated formula for the diffusion 
constant (3.10). Of course, the diffusion constant still varies very smoothly with 7 in 
the physical regime, as shown in Fig. 2. We will confirm the above result by directly 
evaluating the two-point function of the dual current in the next section. 

Given these results and the Einstein relation D = a /x, the susceptibility is easily 
determined to be 



X 



9l 



arctan 



1 + 7 



1/3 



16ttT7 1 /3 V/ -v. LV37 1 ' 3 

Again considering small 7, the susceptibility reduces to 

4ttT / _ 36 



log 



41 



2 7 l/3)c 



X 



2 7 



7 



(3.13) 



(3.14) 
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Figure 2: The charge diffusion constant is plotted versus the coupling 7. The vertical 
dashed lines denote the boundaries of the physical regime, 7 G [—1/12, 1/12] - see discussion 
surrounding eq. 5.14. 



4. Conductivity 

In this section, we calculate the conductivity for the CFT dual to the bulk action (2.6). 
We begin by decomposing the gauge field as 

A io-y 



A a (t, x, y,u)=J ^e^A a (u, q) , (4.1) 



where q • x = — out + q x x + q y y. For convenience and without loss of generality, we 
choose three- momentum vector to be q M = (u,q,0). Further we choose the gauge in 
which A u (u, q) = 0. Then evaluating modified Maxwell's equations (2.7) in the planar 
black hole background (2.4), we find 

g/(3-16nV") 
A + u;(3 + 32u 2 7 f) K ~ ° (4 - 2) 
4n 7 (2f" + uf") L 4 q (3 - 2n 2 7 f) 



All /'(3 - 2m 2 7/") - 2u 7 /(2/" + u/'") „ L 4 w , 4 



,„ f{Z-2u* 1 f")-2u 1 f{2f" + uf") 
/(3 - 2 M 2 7 /") 



L 4 (3c 2 - 3g 2 / - 2^ 2 7 (^ 2 + 2g 2 /) fQ , n r , 
+ P (3 - 2 M 2 7 f ') (45) 
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Now we can use equations (4.2) and (4.3) to decouple equation of motion for A t (u,q): 



A';' + gi (u)A'; + g 2 (u)A' t = 0, (4.6) 

where 

f (9 + 6w 2 7 /" - 64w 4 7 2 f 2 ) + 2n 7 /(15 - 4m 2 7/")(2/" + wf") 
~ /(3-2 M 2 7 f)(3 + 2 M 2 7 f) ' 

^ = r g /»(3 - 2, 2 7 n(3 + I^j ^ + 6 ^ - ^V/") 
+ /(3 - 2u 2 1 f")(-3q 2 L* + 2 U1 {q 2 L% + 4r 2 /')f + 4r V 7 /'/'") 
+8r 2 7 / 2 (3f + 2n 2 7 /" 2 + Quf" + u^f" 2 )) . (4.7) 

At this point, recall that in the analysis of the Maxwell theory in [7], the equations 
of motion for A y (u, q) and A' t (u,q), i.e., the 7 = limit of eqs. (4.5) and (4.6), were 
identical. This was a result of the EM self-duality of this bulk theory. However, clearly 
eqs. (4.5) and (4.6) are no longer identical with nonvanishing 7 , indicating that the 
new interaction in eq. (2.6) breaks the EM self-duality in the present case. We return 
to examine the EM duality in detail in section 6. 

Next we solve eq. (4.5) with an infalling boundary condition at the horizon. Near 
the horizon, we can write A y (u, q) = (1 — u) b F(u, q) where F(u, q) is regular at u = 1. 
Inserting this ansatz in eq. (4.5), we find that b = ±iL 2 uj/(3r ). The ingoing boundary 
condition at the horizon fixes 

L 2 u 

b = -i = —iw, (4.8) 

3r 

where we have defined the dimensionless frequency 

to=-^=. (4.9) 

As we wish to calculate the conductivity with w^O but q = (recall that q is spatial 
momentum along x-direction), we simplify the notation by denoting A a (u,u,q = 0) 
and F(u,ou,q = 0) by A a {u) and F{u). With b given by (4.8), for q = 0, the equation 
of motion for F(u) reduces to 

°~ F + { (1- U 3)(l +4U 3 7 ) -T^) F 

ivo ((1 + u + u 2 )(l + 2u + 2u 2 (3 + 4u + 5u 2 ) 7 ) - i{2 + u)(4 + it + u 2 )(l + 4u 3 7 )t») 
+ (1 - u)(l + u + u 2 ) 2 (l + 4u 3 7 ) 
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To proceed further, we need to recall the relation of the conductivity to the retarded 
Green's function Q yy for the dual current J y : 



a 



-Im 



g w (q) 



(4.11) 



Of course, we wish to calculate Q yy using the AdS/CFT correspondence, following [37]. 
Briefly, integrating by parts in the action (2.6), the bulk contribution vanishes by the 
equations of motion (2.7) and so the result reduces to a surface term. At the asymptotic 
boundary, one has the following contribution for A y 



1 uu 



\ [ d 3 x^g uu gyy(l-8 1 L 2 C uy u y)A y (u,x)d u A y (u^) 
d 3 xA y (u, x)d u A y (u, x) 



2<?4 

2ttT 



m-s-0 



(4.12) 



u->0 



The simple expression in the second line results from explicitly evaluating the expression 
with the black hole metric (2.4) for which 



^uy 



U 



2L 2 



(4.13) 



The Fourier transform of A y is required to compare the above expression with the 
standard AdS/CFT result 



l yy 



d 3 q 1 
(2tt) 3 2 



A y {-q)Q yy {q) Ay(q) 



u->0 



(4.14) 



Hence we can arrive at the usual result, i.e., the coupling 7 makes no explicit appearance 
here, 

4vrT A y (u, -q) d u A y (u, q) 



Gyy(q) 



(4.15) 



3gj A y (u,-q)Ay(u,q) 

Focusing our attention on the case q M = (ou, 0, 0) and adopting the notation introduced 
above eq. (4.10), the retarded Green's function becomes 



Q yy {u,q = 0) 



AnT d u A y (u, to) 



3gj A y (u,u) 



u-+0 



(4.16) 



Then eq. (4.11) yields the conductivity at q = as 



= T~2 lm 



d u A % 



to A, 



(4.17) 



y / 
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Given the above expression, it is straightforward to calculate conductivity for small 
lo analytically and confirm the result (3.12) for cr = cr(w = 0, q — 0) derived in the 
previous section using the membrane paradigm. First, we make a Taylor expansion of 
F(u) in vo and substitute the ansatz F(u) = Fi(u) +voF2(u) into (4.10). Then, we find 
that Fi and F 2 should satisfy the following 

3w 2 (l -4(1 -2w 3 ) 7 ) 



F'i 



F'l 

3m 2 (1-4(1-2m 3 ) 7 ) 



F' 



W 3 )(l 

2i 



4m 3 7) 



-F[ = 0, 



(4.18) 



m 3 )(1 + 4w 3 7) z 1 — u 
i(l + 2u + 4n 2 7(3 + 4u 



+ 



5m 2 )) 



F 1 = 



(4.19) 



(1 - u 3 )(l + 4m 3 7 ) 

After solving eq. (4.18) for F\, we can fix one of the integration constants demanding 
that Fi is regular at the horizon. This yields Fi(u) = C, where C is an arbitrary 
constant. Given Fi, we solve eq. (4.19) for F 2 . In this case, we fix the two integration 
constants by imposing the following two conditions: First, F 2 is regular at the horizon. 
Second, we normalize F(u) such that its value at the horizon is independent of vo, i.e., 
F 2 (u — 1) = 0. The final result is given by 



F 2 (u) 




arctan 



2v / 3 2 1/3 7 2/3 ( arctan 



2u 



V3 J 



log 



1 + U + 7T 



(4.20) 



2 2 2 / 3 m 7 1 / 3 

~7! 



arctan 



1 _ 2 2 2 / 3 7 1 /3 



+ 2 1 / 3 7 2 / 3 log 



(i + 2 2 /y/ 3 ) 
1 + 47 



2 1 / 3 7 2 / 3 log 



(l + n2 2 /y/ 3 ) 3 
1 + 4u 3 7 



Now we can simply use A y (u) ~ (1 — u) b (Fi(u) + voF 2 (u)) in eq. (4.17), take the limit 
to — > and find 



CTO = -2(1 

91 



4 7 ) 



(4.21) 



which agrees with our previous result (3.7). 

To study frequency dependant conductivity, we must solve eq. (4.10) numerically. 
Our numerical integrations run outward from the horizon and so we need to fix the 
initial conditions at u = 1. To determine the latter we solve eq. (4.10) for « <C 1, 
finding 

ivo(i + 2t» + 87(2^ + t»)) 



F(u) 



u 



(4.22) 



(l + 4 7 )(i + t») 

Numerical integration is used to determine F(u) out to the boundary at u = for fixed 
values of vo (and 7) and then we use the complete solution A y {u) = (1 — u) b F{u) and 
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eq. (4.17) to calculate conductivity <r(t»). In figure 1, we show our results for various 
values of coupling constant 7. 



5. Bounds on the Coupling 

In this section, we find the constraints that are imposed on the coupling 7 by demanding 
that the dual CFT respects causality, following the analysis described in [31, 29, 35]. 
We also examine if there are any unstable modes of the vector field, as discussed in 
[38, 30], which would result in our calculations of the charge transport properties being 
unreliable. From a dual perspective, such unstable modes indicate that the uniform 
neutral plasma is an unstable configuration in the dual CFT. 

To examine causality, the first step is to re-express the equations of motion of the 
two independent vector modes, i.e., eqs. (4.5) and (4.6), in the form of the Schrodinger 
equation. We begin by considering eq. (4.6). Recall that we are working in the gauge 
where A u (u, q) = and we have chosen q M = (u,q,0). Now if we make a coordinate 
transformation to z(u) such that 

1 — u 

and write A' t (u, q) = G\(u) ipi(u, q) where 

n>(\ 6n 2 7 (5 + 8n 3 7 ) 

G, M r— : — — Gl(Ul=(J, (5.2) 

1V ; 1 - 4n 3 7 (l - 64n 3 7 ) u ' v ' 

then eq. (4.6) takes the form 

-d^z) + V(z)if>i(z) = vo 2 i>(z) . (5.3) 

In this Schrodinger form, the effective potential V(z) can be expressed in terms of u as 

V( U )=c f 2 V (u) + V 1 (u) J (5.4) 



where 



(l-n 3 )(l + 4n 3 7 ) R . 

V ° {U) = (1 - 8n 3 7) (5 - 6) 

2 M (l- M 3 )7(2-2 M 6 7-5 W 3 (l + 27)) 

Vl{U) ~ 3(l + 4n 3 7 ) 2 • { ] 
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It is easiest to consider the limit cj — > oo, in which case one can solve for ipi in a 
WKB approximation [31]. In this limit, Vq(u) will dominate the potential and we want 
to examine how its properties change as 7 is varied, e.g., following [29]. In Fig. 3, we 
have plotted potential Vq(u) for various values of 7. We observe that if 7 is too large, 
the potential develops a maximum with Vo tmax > 1 at some point between u = and 
u — 1. In that case, there will be 'super-luminal' modes with t»/q = uj / 'q > 1 indicating 
that causality is violated in the dual CFT [31, 29]. One can easily verify that this 
new maximum appears for 7 > 1/12 by examining the behaviour of Vq(u) near the 
boundary, i.e., near u = 0, where eq. (5.6) yields 

V {u) ~ 1 - (1 - 12 7 ) m 3 + ••• . (5.8) 




Figure 3: Left: Vq(u) for longitudinal At mode for various values of 7. Right: Wq(u) for 
transverse A y mode for various values of 7. Consistency conditions discussed in the text 
for the longitudinal mode imply 7 £ [—1/4,1/12] in Vq. Similarly for transverse mode, 
7 G [-1/12,1/8] in Wq. 



Next we turn to the transverse vector mode satisfying eq. (4.5). As above, we make 
a change of coordinate to z(u) satisfying eq. (5.1) and we write A y (u) = G2(u) ip2{u) 
where 

G>(u) + Y ^-G 2 (u)=0. (5.9) 
With these choices, eq. (4.5) reduces to the desired Schrodinger form 



-d 2 Mz) + W{z)^ 2 {z) = xx> 2 ^ 2 {z) , 



(5.10) 



where 



W(u) 
W (u) 

Wx(u) 



<CW (u) + Wi{u) with 
(1 -m 3 )(1 -8w 3 7 ) 

(1 + 4u 3 7 ) ' 
2u{l - w 3 ) 7 (2 - 5m 3 + 2 7 w 3 (l - 7u 3 )) 
3(1 + 4w 3 7 ) 2 



(5.11) 
(5.12) 

(5.13) 



We again consider the WKB limit where Wq dominates the potential. The shape of this 
potential is also shown in Fig. 3 for various values of 7. Examining the potential (5.12) 
as above, we find that a maximum develops for 7 < —1/12, indicating that causality 
is violated in the dual CFT in this regime. 

Combining the results from both modes, we find that the dual CFT is only consis- 
tent (i.e., respects causality) if 



1 1 

— < 7 < — 
12 ~ ' ~ 12 



(5.14) 



We also note that these bounds coming from the violation of micro-causality precisely 
match the bounds derived for the dual parameter in the CFT derived in [24, 34]. 
There, various thought experiments were proposed to constrain CFT's in four dimen- 
sions. However, their discussion is readily adapted to the three dimensions, as we 
consider here. The relevant experiment consists of first producing a disturbance, which 
is localized and injects a fixed energy, with an insertion of the current e l Ji, where e l is 
a constant (spatial) polarization tensor. Then one measures the energy flux escaping 
to null infinity in the direction indicated by a unit vector n: 



Sin) 



"+00 

lim r ' 

r— >+oo 



dt T*i(i,rn) n 1 . 



(5.15) 



The final result takes the form 

(0|(£*-jt)£(n) (s-J)\0) 



<£(n)> 



(0\(e* ■ Jt) (e ■ j)\0) 



E 
2^ 



1 + 0,2 



£ • n 



E 
2^ 



1 + a 2 ( cos 2 9 



(5.16) 



where E is the total energy and 9 is the angle between the direction n and the polar- 
ization £. The structure of this expression is completely dictated by the symmetry of 
the construction and the (constant) coefficient a 2 is a parameter which characterizes 
the underlying CFT. Given eq. (5.16), it is clear that a 2 is related to the parameters 
appearing in the general three-point correlator {T ab (x)J c (y)Jd(z)) - see discussion in 
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section 7. Now, the interesting observation of [24] was that if the coefficient a 2 becomes 
too large, the energy flux measured in various directions will become negative. Hence 
demanding that the energy flux should be positive in all directions for a consistent CFT 
leads to the constraints 

-2<a 2 <2. (5.17) 

Of course, to relate this result to that in eq. (5.14), we must find the relation between 
a 2 for our holographic CFT and the bulk coupling 7. The simplest approach is to use 
the AdS/CFT correspondence to examine the bulk dual of the thought experiment pre- 
sented above. As noted above, in calculating the flux expectation value in eq. (5.16), 
we are essentially determining a specific component of the three-point function of the 
stress tensor with two currents. Hence in our holographic description, we must intro- 
duce an appropriate metric fluctuation h^ u and two gauge field perturbations in the 
AdS4 bulk, which couple to the boundary insertions of T a b and J a . We then evaluate the 
on-shell contribution for these three insertions with the action (2.6). We do not present 
the details here, as the analogous calculations for d = 4 are presented in Appendix D 
of [24] - the interested reader may also find the discussion in the first reference in [25] 
useful. In the end, the holographic calculations yield a very simple final result 

a 2 = -24 7 (5.18) 

and hence we find the bounds in eqs. (5.14) and (5.17) are equivalent. 

Next we turn to possible instabilities in the neutral plasma. If we examine the 
potential Vq in more detail, we find that another interesting feature develops for 7 < 
— 1/4. That is, the potential develops a minimum at some radius close to the horizon 
where Vq(u) < 0. The appearance of this potential well can be verified analytically by 
expanding Vq{u) near u — 1, 

V (u) ~ 3^il(l - u) + ■■ ■ . (5.19) 
1 — 07 

While Vq always vanishes at u — 1, we see that for 7 <— 1/4, Vo<0 immediately in 
front of the horizon indicating the presence of the negative potential well there. In the 
WKB limit, this potential well leads to bound states with a negative (effective) energy, 
which correspond to unstable quasinormal modes in the bulk theory [38]. While these 
modes do not signal a fundamental pathology with the dual CFT, they do indicate 
that the uniform neutral plasma is unstable in this regime. Hence our calculation of 
the conductivity would be unreliable here. Of course, our previous constraints (5.14) 
have already ruled out 7 < —1/4 as being physically interesting and so we need not 
worry about these instabilities. 
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On the other hand, one may worry that additional instabilities will appear outside 
of the WKB regime, considered above. In particular for small momentum, the effective 
potential will also receive an important contribution from V\(u). We find that for 
7 G (—1/4,0), V\{u) also develops a negative minimum close to the horizon and so 
there might be some unstable modes in the plasma in this regime as well. We have 
plotted the potential Vi(u) for various values of the coupling constant 7 in Fig. 4. 
While the WKB approximation may be less reliable in this regime, the analysis in [38] 
suggests that it is sufficient to determine the appearance of unstable modes. According 
to WKB approximation, a zero energy bound state can appear in this potential well 
for 



n 



7T 



20 



dzyJ-V^z) 

3du 



-V x {u) 



(5.20) 



where n is a positive integer and the integration is over the values of u for which the 
potential is negative. A plot oin = I/n + 1/2 is given in Fig. 4. We see that h reaches a 
maximum value of approximately 0.86, implying that the potential well is never able to 
support a negative energy bound state. Hence we conclude that there are no unstable 
modes in this low momentum regime. 




-0.20 -0.15 -0.10 -0.05 



Figure 4: Left: V\{u) for various values of 7. Right: n = I/ir + 1/2 plotted versus 7 - see 
eq.(5.20). Here potentials Vi(u) is plotted for various values of 7. We see that a negative dip 
appears in V\{u) close to the horizon for 7 S (—1/4, 0). We have also plotted n in this range 
of 7 and the plot clearly indicates that it always remains less than one. 



While we have discussed both small and large momenta limit of our effective poten- 
tial V(u), one may still imagine that instabilities can still arise at some finite momenta. 
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However, such a possibility can be eliminated by considering the structure of our com- 
plete potential V(u). That is, for any finite momenta and for 7 G [—1/4, 0], the negative 
dip in potential V(u) is smaller than the dip in V\{u) because of the positive contribu- 
tion coming from Vq(u). Hence there are no instabilities coming from the longitudinal 
vector mode in the regime (5.14) of physical interest. 

Of course, one must also consider possible instabilities in the transverse vector 
mode. In this case, examining the potential Wo, we find that a negative minimum 
again develops for 7 > 1/8. So again instabilities appear in the large momentum limit 
but only for values of the coupling outside of the physical regime (5.14). As above, 
one can also consider the low and finite momentum regimes, however, again one finds 
that there are no additional instabilities in the physical regime. Hence although both 
the transverse and longitudinal modes of the vector exhibit instabilities, these only 
appear in a regime where our previous constraints already indicate that the CFT is 
pathological. 

Examining eq. (5.19), one sees that the potential Vq is also negative in front of 
the horizon for 7 > 1/8 (as well as for 7 < —1/4, as discussed above). However, 
this behaviour is not indicative of a negative potential well in this case. Rather a 
closer examination of the full potential (5.6) shows that a simple pole appears at u — 
l/(27 1 / 3 ), which lies in the physical interval < u < 1 for 7 > 1/8. The potential Wq 
exhibits a similar behaviour for 7 < — 1/4. The analysis and physical interpretation of 
the modes in this case are more elaborate along the lines of that given in [25] . However, 
we do not consider these issues further here since our previous constraints (5.14) already 
indicate that 7 > 1/8 and 7 < —1/4 are outside of the physically viable regime. 

6. EM Self-Duality Lost 

In this section, we examine in more detail the loss of electromagnetic (EM) self-duality 
for the U(l) gauge theory defined by the bulk action (2.6). Recall from [7] that this 
EM self-duality was the key property of the standard four- dimensional Maxwell theory 
which lead to the simple relation: 

K T (u), q) K L (co, q) = constant , (6.1) 

where K T and K L are the scalar functions determining the transverse and longitudinal 
components of the retarded current-current correlator - see appendix A for further 
discussion. As a result, the conductivity (at zero momentum) was a fixed constant for 
all values of uj/T. In examining the explicit equations of motion, (4.5) and (4.6), we 
already noted that self-duality is lost in the new theory. However, in the context of 
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any U{1) gauge theory, one can think of EM duality as simply a change of variables in 
the corresponding path integral. Even if our new gauge theory (2.6) is not self-dual, 
we can still implement this change of variables and construct the EM dual theory, as 
we will demonstrate below. 

We begin by introducing a (vector) Lagrange multiplier B a in the generalized action 
(3.1) as follows 

/ = J d 4 xV^9 (-^F ab X abcd F cd + l -e abcd B a d b F cd ^ . (6.2) 

Here e abcd is totally antisymmetric tensor, with £0123 — \f~9- The fundamental fields 
in the path integral for this action are the two-form F ab and the one-form B a . Now the 
EM duality comes from simply treating the integration over these fields in two different 
orders. 

If we evaluate the path integral by first integrating over the Lagrange multiplier 
B a , the latter integration enforces the Bianchi identity on the two-form F a b, i.e., 

e abcd d b F cd = 0. (6.3) 

If F ab is to satisfy this constraint, 3 then on a topologically trivial background, it must 
take the form F ab = d a A b — d b A a . Hence the remaining path integral reduces to the 
'standard' gauge theory where the fundamental field is the Maxwell potential A a with 
generalized action given in eq. (3.1). 

Alternatively, one can perform the path integral over the two- form F ab first. In this 
case, we first integrate by parts in the second term in the action (6.2) 

1 = J d A x^g (y-l^F ab X ahcd F cd + \e abcd F ab G c )j . (6.4) 

where we have defined the new field strength G ab = d a B b — d b B a . We can now shift the 
original two-form field to 

F ab = F ab - 9 l{X-') abcd e^G ef (6.5) 

where X~ l is defined by 

^{X-^X^ = 1^ . (6.6) 



3 Note that we are justified in using ordinary (rather than covariant) derivatives both here and in 
the action (6.2) because of the antisymmetry of the indices. 



- 20 - 



Recall the definition of I ab cd given in eq. (3.3). With this shift, one has a trivial Gaussian 
integral over the field F ab after which one is left with the path integral over the one-form 
B a with the action 

1 = J d 4 xV=9 (^-^X abcd G ab G c )j , (6.7) 

where g\ = 1/gf and 

V cd _ ef I v — IN qh _ cd 

A a fe — ~^ £ ab (A )ef £gh 

= (X-\ b - + \(X-X f "l^ (6.8) 

\ I V — 1\ cec d I \r — 1\ dec c ( V — 1\ cer cZ , / v — 1\ dec c"| 
- [(A J ae Of, - (A J ae d b - (A Jfce d a + (A Jfe d a J . 

In the second equality above, the two e-tensors have been eliminated with the four- 
dimensional identity of the form e a bcd^ 9h = —(S a e S^S c 9 5d h + • • • ) - note that the 
overall minus sign appears because we are working in Minkowski signature. Hence B a 
now plays the role of the gauge potential in the EM dual theory with the action (6.7). 

The relation between the gauge fields in the two dual EM theories is implicit in the 
equations of motion for F a b or F a &. From eq. (6.5), we see that setting F ao = yields 



,2 



* ab 



fh ( V -l\ cd 

4 



(X- 1 ) "ej'G*. (6.9) 



Recall that in the usual Maxwell theory, X takes the simple form given in eq. (3.3). 
In this case, X -1 = X and one can easily show that eq. (6.8) also yields X ab cd = I a b cd - 
Hence for the Maxwell theory, the form of the two actions, (3.1) and (6.7), as well 
as the corresponding equations of motion for A a and B a , are identical. This is then 
a demonstration that the Maxwell theory is self-dual. Further, the duality relation 
between the two field strengths in eq. (6.9) corresponds to the usual Hodge duality, as 
expected for this case. 

Of course, in general, we will find that X ^ X and so this self-duality property is 
lost. That is, the form of the action and the equations of motion in the original theory 
and its dual now have different forms, i.e., 

V b {X abcd F cd ) = and V b (x abcd G cd ) = . (6.10) 

For the action of interest (2.6), X is given in eq. (3.4) and at least in a regime where 
we treat 7 as small, we can write 

(X- l ) ab cd = I ab cd + ^L 2 C ab cd + 0( 7 2 ) • (6.11) 
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Further because of the traceless property of the Weyl tensor, one finds 



X ab cd = (X- 1 ) afo cd + 0( 7 



(6.12) 



With the change in sign of the order 7 contribution between eqs. (3.4) and (6.11), it is 
clear that our gauge theory is no longer self-dual. 

Actually given the planar black hole background (2.4), it is straightforward to 
calculate X~ x exactly. First, we define a six-dimensional space of (antisymmetric) 
index pairs with, i.e., A, B £ {tx,ty,tu,xy,xu,yu} - note both the ordering of both 
the indices and the index pairs presented here. Then X given in eq. (3.4) becomes a 
diagonal six-by-six matrix 



X, 



diag (1 + a, 1 + a, 1 



2a, 1 



2a, 1 + a, 1 



a 



(6.13) 



where a = 47 u 3 . Since X is a diagonal matrix, X^ 1 is also a diagonal matrix whose 
entries are simply the inverses of those given in eq. (6.13). Note that a takes its 
maximum value at the horizon u = 1, i.e., a max = 47. Hence we must constrain 
-| < 7 < I in order for the inverse to exist everywhere in the region outside of the 
horizon. Of course, it is not a coincidence that the effective Schrodinger equation in 
section 5 became problematic (i.e., the effective potentials contained a pole) precisely 
outside of the same interval. In any event, the physical regime (5.14) for 7 determined 
in section 5 lies well within this range. 

Using this notation and the background metric (2.4), e a t, cd becomes the following 
'anti-diagonal' six-by-six matrix 



B 



rpf 
' L 2 



rpf 
L 2 



ro 



rof 



rof 



ro_ 
L 2 



(6.14) 



Combining these expressions, we can easily evaluate the duality transformation (6.9), 
which is expressed using the new notation as 



Fa = 9 1 [X-') A B e B c G 



c 



(6.15) 
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The final result is 



tx ~ 1 + a L 2 u ' ty ~ 1 + a L 2 u ' 



tu 



r ^-G, s/ , = -—^-^G^, ( 6 - 16 ) 

2 r2 2 r2 

3:11 l + «r„/ ty ' yu l + ar f tx - 

This duality transformation gives us a precise analytic relation between the original 
gauge field A a and that, B a , in the EM dual theory. Of course, it would be less 
straightforward to express these duality relations in a covariant construction using the 
Weyl curvature tensor. 

As discussed in [7], from the perspective of the boundary field theory, we can 
describe the CFT in terms of the original conserved current J a (dual to the bulk vector 
A a ) or a new current J a (dual to B a ). In the case of the Maxwell theory, the EM self- 
duality means that both currents have identical correlators. In the present case, where 
EM self-duality is lost, the correlators still have a simple relation which is summarized 
by 

K T (cu,q) K L (u,q) = 1, (6.17) 
K T (u,q) K L (u,q) = 1. 

The detailed derivation for these relations can be found in appendix A. The self-dual 
version of eq. (6.17), with K = K, appeared in [7]. However, the conventions for the 
EM duality transformation were different there, i.e., they chose §4 = g^. This choice 
changes the normalization of the dual currents and so changes the constant on the right- 
hand-side of eq. (6.17) to (<?4)~ 4 . In any event, these relations imply, the longitudinal 
correlator in one theory is traded for the transverse correlator in the dual theory, as 
reflected in eq. (6.16). Notably, eq. (6.17) has precisely the same form as that obtained 
from general considerations of particle-vortex duality, but without self-duality, in the 
condensed matter context, as we review in the following subsection. 

6.1 Particle- Vortex Duality 

Above, we discussed EM duality as a change of variables which allows us to formulate 
the bulk theory in terms of two different gauge potentials. This reformulation of the 
bulk theory implies that the boundary CFT can be developed in terms of two 'dual' sets 
of currents, whose correlators are simply related using eq. (6.17). As noted in [7], the 
latter is reminiscent of the structure of the correlators in systems exhibiting particle- 
vortex duality. The discussion there focused on self-dual examples, however, the latter 
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is an inessential feature to produce eq. (6.17), as we illustrate with the following simple 
example - see also appendix B of [7]. 

Consider the field theory of a complex scalar z coupled to a U(l) gauge field 



S 



dxdt 



-12 , „.|„|4 , 1 



1(8^ — iA„) z\ + s\z\ + u\z 



2e 2 



{e^ x d u A x )' 



(6.18) 



We now look at the structure of the conserved U(l) currents of S, and their correlators. 
For simplicity, we will restrict our discussion to T = to make the main point in the 
simplest context. There is a natural generalization to T > 0, which is needed to obtain 
the full structure of the relationship in eq. (6.17), and which was discussed in [7]. 
The theory S has the obvious conserved U(l) current 



J, = -z* [(d, - iA,)z\ - - [(d, - iA^z*} z 



(6.19) 



Because of current conservation, we can write the two-point correlator of this current 
in the form (reminder, we are at T = 0) 



'Mp)Ju(-p)) 



fJbV 



Pjjh 

p2 



(6.20) 



Here, we note that this correlator has been defined to be irreducible with respect to 
the propagator of the photon, A^. 

The theory S has a second conserved U(l) current; this is the 'topological' current 



^ ~ 27T 



d v A> 



(6.21) 



We can interpret as the current of dual set of particles which are the Abrikosov- 
Nielsen-Olesen vortices of the Abelian-Higgs model in eq (6.18). Each such vortex 
carries total A^ flux of 2n, and hence the prefactor above. Indeed, there is a dual 
formulation of the theory in eq (6.18) in which the vortices become the fundamental 
complex scalar field 'z: 



S 



dxdt 



|<9 M £| 2 + s~|^| 2 + u\z\ 



(6.22) 



This dual theory has no gauge field because the vortices of S only have short-range 
interactions. The particle number current of this dual theory is the same as that in 
eq. (6.21) 

Jfj, = -z'*d ll z — rdftZ* z. (6.23) 

% % 
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Now, returning to the perspective of the original theory S in eq (6.18) and the U(l) 
current in eq. (6.21), we can write the two-point correlator of in the general form 

(W(-ri)-^(^-^) p2/ /_^ 2) (^4) 

where XQo 2 ) is the photon self-energy. 

The photon couples linearly to the current J^, and so the photon self energy is 
clearly the irreducible J M correlator, and so 

S(p 2 ) = ^K(p 2 ) (6.25) 

Also as p 2 — > in IR, we have S(p 2 ) ^> p 2 /e 2 - recall that here we are assuming the 
spacetime dimension d = 3. So we have 

(j»(p)U-P)) - - - Pj ^A v9K(P 2 ) (6-26) 
where from eqs. (6.24) and (6.25) 

K(p 2 )K(p 2 ) = ^. (6.27) 

This result is clearly the T = analog of eq. (6.17). It is easily generalized to T > 
0, after separation into transverse and longitudinal components, but we refrain from 
presenting those details here. 



7. Discussion 



Our main results for the frequency dependence of the conductivity without self-duality 
were given in Fig. 1, and we presented a physical interpretation in Section 1. For 7 > 0, 
the results had a qualitative similarity to that expected from a Boltzmann transport 
theory of interacting particles, while for 7 < the results resembled the Boltzmann 
transport of vortices. 

We will now discuss other aspects of these results. We also see from Fig. 1 that 
the large frequency limit is unaffected by the new coupling, i.e., a{u = 00) = 1/gf. 
We can understand this result from the fact that the Weyl curvature vanishes in the 
asymptotic region of the black hole region and so the new interaction in eq. (2.6) has 
no effect there. 

Further, we have 

<t(u = 0) . , 

-y ^ = l + 4 7 7.1 

a{oo = 00) 
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and so this ratio varies between 4/3 and 2/3 in the allowed physical regime given in 
eq. (5.14). Thus the allowed range of variation in the conductivity by non-self-duality 
is smaller than 33% and can have either sign, in our model. This should be contrasted 
from the large variation obtained from the weak-coupling Boltzmann analyses. In the 
e = 4 — d expansion (where d is the spacetime dimension), it was found that generically 
[6] 

a{uj = oo) e z 

Similarly, in the large N expansion (where N is the number of components of a vector 
(and not matrix) field), we have [8] 

^\~JV. (7.3) 
a(uj = oo) 

In both cases, the ratio becomes large in the regime of applicability of the analysis. 
Thus the AdS/CFT analysis gives a useful result for this ratio in the complementary 
limit of very strong interactions. 

Also note that the conductivity in Fig. 1 does not vary monotonically, rather it 
seems there is an extremum at u> ~ 2irT. For 7 > 0, this oscillation is as anticipated 
from Drude-like considerations of particle transport in ref. [6] , and for 7 < we argued 
in Section 1 that such an oscillation is obtained from Drude-like vortex transport. 
Recall that in the AdS/CFT correspondence, particle- vortex duality in the boundary 
theory is realized as EM duality in the bulk theory. Hence we can make the previous 
point explicit for our holographic model using the formalism developed in Section 6. 
That is, for any given value of 7, we can explicitly construct the EM dual theory and 
evaluate the conductivity. In Fig. 5, we have plotted the resulting conductivities for 
the original bulk theory and the EM dual theory for 7 = ±1/12. As expected, for 
7 = —1/12, the conductivity of the dual theory exhibits a Drude-like peak at small 
u. For 7 = 1/12, a similar peak appears for the original theory while the EM dual 
theory exhibits a dip in the conductivity at small u. For either value of 7, the figure 
also illustrates that the conductivities of the two dual theories are not precise inverses 
of one another, except for u — > 0, 00. This occurs because the function K T (u,q = 0) 
is only precisely real in the latter limits. 

Oscillations in the conductivity similar to those in Fig. 1 were observed in [39]. 
The latter studied the transport properties of currents on a three-dimensional defect 
immersed in the thermal path of a four-dimensional superconformal gauge theory. The 
holographic bulk theory consisted of probe D-branes embedded in AdSs x S 5 and the 
oscillations were an effect of stringy corrections to the usual D-brane action. Implicitly, 
the four-derivative interaction considered there would have been a linear combination of 
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0.5 - 



Figure 5: The (dimensionless) conductivity a = g\a is plotted versus the (dimensionless) 
frequency w = uj/{AttT) for various values of 7. The solid curves correspond to the same 
conductivities displayed previously in Fig. 1 — red for 7 = 1/12 and blue for 7 = —1/12. 
The dashed curves show the conductivity calculated from the EM dual theory for the same 
values of 7. 



the 05,6,7 terms in eq. (2.5). In this previous setting, the calculations were perturbative 
and the oscillatory contribution to the conductivity was suppressed by a factor of A -1 / 2 
relative to the constant term produced by the Maxwell action on the brane - as usual, 
A denotes the 't Hooft coupling of the four-dimensional gauge theory. 

Section 2 provided some motivation for introducing the new four-derivative inter- 
action in eq. (2.6). However, there was a certain liberty in choosing the precise form 
of the curvature in this interaction. From a certain perspective, the following vector 
action may be preferred: 



IL 



d 4 Xy/—g 



F ab F ab + aL 2 



[R 



Fab rpcd 



AR 



ab 



F ac F b c + RF ab F, 



ab] 



(7.4) 

The advantage of the higher-derivative term above is that it produces second-order 
equations of motion for both the gauge field and metric in any general background. We 
can think of this term arising from Kaluza-Klein reduction of Gauss-Bonnet gravity in 
five-dimensional space-time [40]. Now the generalized Maxwell's equations are 



V a [F ab -AaL 2 {R ab cd F cd 



2R ac F c " + 2R bc F c a + RF 



Tiab 



0. 



(7.5) 



Before considering the charge transport for this theory, we note that AdS vac- 
uum and the neutral black hole (2.2) remain unmodified with this choice of the four- 
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derivative interaction. In particular then, for the black-hole background, we still satisfy 
the vacuum Einstein equations, i.e., R a b = —3/L 2 g a b. Further, the Reimann curvature 
tensor R a bcd is related to the Weyl tensor C a bcd by 

Rabcd = C a bcd + 9a[cRd]b ~ 9b[cRd\a ~ ^Rga[cgd]b ■ (7.6) 

By substituting these relations into eq. (7.4), we find that the action becomes 

1 + 8a f 4 . _/ 1 b a 2 T^abrpcd 



C = J d xV= ~ 9 {-^ F + T+Sa ) ■ < 77 > 

Hence, this expression for action is identical to eq. (2.6) if we identify the couplings: 

2 _94_ and (7. 8 ) 
y4 l + 8a ' l + 8a v ; 

Hence in the neutral plasma, all of the charge transport properties of the new theory 
are identical to those found in the main text, as long as we make this identification of 
the couplings in the bulk gauge theory. For example, we have explicitly applied the 
analysis of section 5 to the new action (7.4) and found this produces the constraints 
—1/20 < a < 1/4. One can easily verify that this range precisely matches that in 
eq. (5.14) for 7 using the identification of the gauge theory couplings in eq. (7.8). 

It would be interesting to examine charged black holes in this new theory (7.4). 
Beyond analyzing the effects of adding a chemical potential in the boundary CFT, it 
would be interesting to examine the so-called "entropy problem" in this theory. That 
is, at zero temperature, charged black holes still have a finite horizon area for the 
Einstein-Maxwell theory in the bulk and hence the dual CFT has a large entropy even 
at T = but nonvanishing chemical potential. It would be interesting to determine 
how this feature found in simple holographic CFT's is affected by the introduction of 
the new higher derivative bulk interaction in eqs. (2.6) or (7.4). Such investigations 
would require numerical work that would be greatly facilitated by having second-order 
equations, as produced by the above action (7.4). 

As discussed in the introduction, we are following a program of expanding the uni- 
versality class of the holographic CFT by introducing new higher-derivative interactions 
to the bulk action. The simplest way to characterize the effect of the new interactions 
is to examine the changes which are produced in the vacuum n-point functions in the 
CFT. As alluded to above, the Weyl curvature vanishes in AdS space and so we may 
infer that in the vacuum of the dual CFT (with vanishing temperature and charge 
density), there are no changes to any of two-point functions, i.e., {J a (x)Jb(y))o and 
{T a b(x)T cd (y)) , where the subscript indicates the two-point functions are evaluated 
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in the vacuum or at T = 0. That is, the two-point functions are independent of 7 in the 
vacuum. In particular then, the charge transport properties of the holographic CFT 
must be independent of 7 at T = 0. On the other hand, recall the simple 7 dependence 
which appears in eq. (3.12) for the conductivity at u — 0. Clearly, this means that the 
limits, T — > and u — > 0, do not commute, as was also emphasized in ref. [6]. 

As described in [24, 34], the key effect of the new bulk interaction (2.6) is to 
modify the three-point correlator {T a b(x)J c (y)Jd(z)) Q . One can show that in any CFT, 
conformal symmetry will completely fix this three-point function between the stress 
tensor and two conserved currents up to two constant parameters [26]. One of these 
parameters vanishes in the holographic dual of an Einstein-Maxwell theory. However, 
this extra parameter is nonvanishing for the CFT dual for our extended theory with 
7^0. In particular, as discussed in section 5, the parameter ai in eq. (5.16) is only 
nonvanishing in the boundary CFT when 7 7^ 0. Of course in a thermal bath, the 
expectation value of the stress tensor is nonvanishing. Hence it should be possible to 
use the previous three-point function to infer the leading 7 modification to the two- 
point correlator {J a Jb)T a t finite temperature, e.g., with an approach similar to that 
considered in [41]. In principle then, such a (perturbative) calculation in the CFT 
should already indicate that self-duality is lost. 

Above, we discussed the behavior of the conductivity, which is related to the current 
correlator at zero momentum. We also studied the full momentum dependence of these 
correlators and obtained the duality relation in eq. (6.17), which applied in the general 
case without se//-duality. Remarkably, this has the same form as that obtained by 
applying particle- vortex duality to a (2+l)-dimensional field theory of a single complex 
scalar, as we reviewed in section 6.1: note that this theory is not self-dual (and self- 
duality is not expected in general, except for a particular theory with two complex 
scalar fields [7, 42]). In the single scalar field discussed in section 6.1, K > 

characterize the transverse/longitudinal components of the two-point correlations of the 
current of the scalar particles - see eq. (A. 4) - while K T,L characterize the corresponding 
quantities of the vortex current. 

Of course, the constants on the right-hand-side of eqs. (6.17) and (6.27) are seen to 
be different. In both cases, this constant depends on the conventions used to normalize 
the currents and a new normalization would change the constant in either model. Hence 
one may ask if these relations can be expressed in a way which removes this ambiguity. 
As we will show, one possibility is to replace eq. (6.17) by 

K T (io,q)K L (u J ,q) = aoa 0l (7.9) 
K T (u,q) K L (u,q) = a a , 

where ctq is the conductivity at zero momentum and zero frequency and <5"o is the same 
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quantity for the dual currents. For our holographic model, gq was given in eq. (3.12) 
and given the discussion in section 6, it is a simple exercise to show that <Jq&q = 1. 
Hence, in this case, we easily recover eq. (6.17) from eq. (7.9) above. However, the latter 
equation applies quite generally as we will now show: First, given the expressions for 
a and K T in eqs. (4.11) and (A. 6), respectively, it is straightforward to show that 
lim^o K T (u,q = 0) = cr . Further with vanishing momentum, K L (u,q = 0) = 
K T (u, q = 0) and hence we also have lim^o K T (u, q = 0) = cr . Now if we know that 
the product K T K L is constant, we can evaluate the constant at vanishing momentum 
and vanishing frequency and then our discussion leads us to write eq. (7.9). This 
expression will apply independent of the conventions used to normalize the currents 
and applies equally well for the field theory examples considered in section 6.1 and in 
[7] as for our holographic model. 

The holographic relation of EM duality in the bulk and particle-vortex duality in 
the boundary theory was first noted in [7, 43] and the effect of this bulk transformation 
on the boundary transport properties was further studied in [44] - see also [39, 45]. 
Particle-vortex duality can be extended to an SL(2, Z) action on three-dimensional 
CFT's [43, 46] and the holographic realization of these group transformations on the 
bulk theory was discussed in [43]. In particular, the S transformation corresponds to 
applying EM duality in the bulk. To discuss the T transformation, the bulk action 
must be extended to include a #-term and acting with the T generator corresponds 
to making a 2n shift of 9. Of course, implicitly or explicitly, the previous holographic 
discussions assumed a standard Maxwell action for the bulk vector. It would be in- 
teresting to extend this discussion of the full SL(2, Z) action to the generalized action 
(3.1) introduced in section 3. Associating the S generator with EM duality as in [43], 
one can easily verify that S 2 = — 1 using eq. (6.9). To include the T generator, we 
would need generalize X to include parity violating terms, i.e., nonvanishing Z2 {u) and 
zs{u) in eq. (A. 23). We leave this as an interesting open question. 

To close, we wish to emphasize that our investigation here has considered a simple 
toy model and one should be circumspect in interpreting the results of our analysis. 
While string theory will generate the higher derivative interactions in our action (2.6), 
it certainly also produces many other higher order terms which schematically take the 
form R n F 2 . For example, some such terms, were explicitly constructed (amongst many 
others) and studied in [47]. Any terms with this schematic form would still fall in the 
class of our general action (3.1) and so modify the charge transport properties in a 
similar way. A key feature of our model was that we were able to identify physical 
restrictions which constrained the new coupling 7 to fall in relatively narrow range 
(5.14). As a result, the conductivity remained relatively close to the self-dual value. 
Our expectation is that similar restrictions appear for general string models, however, 
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finding more comprehensive physical constraints in this context remains an interesting 
open question [48]. As seen here and elsewhere [24, 29, 30, 34], the interplay between 
the boundary and bulk theories in the AdS / CFT correspondence is beginning to provide 
new insights into this question. 
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A. Retarded Green's functions and EM duality 

In this appendix, we find the retarded Green's functions of currents in the boundary 
field theory for finite frequency and finite momentum and further we examine the 
relationship between the Green's functions in the two theories related by EM duality 
in the bulk. In this discussion, we work with the general vector action (3.1) and its EM 
dual (6.7). Recall that the relation between coefficients X and X appearing in these 
two actions is given in eq. (6.8) and the field strengths in the two theories are given by 
F a b = d a Ab — dbA a and G a b = d a Bb — dbB a , respectively. Further the duality relation 
between these two field strengths is given in eq. (6.9). 

For simplicity, we will begin by assuming that X a b cd is diagonal in the six-dimensional 
space defined by the antisymmteric index pairs 

A, B e {tx,ty,tu,xy,xu,yu} . (A.l) 

This property holds for the specific theory (2.6) studied in the main text, as shown in 
eq. (6.13). We comment on more general cases at the end of the appendix. Given this 
assumption, we write 

X A B = disig(X 1 (u),X 2 (u),X 3 (u),X A (u),X 5 (u),X 6 (u)) . (A.2) 

Further rotational symmetry in the xy-plane would restrict this ansatz with X\{u) = 
X 2 (u) and X 5 (u) = X e (u). However, we leave this symmetry as implicit, since it is not 
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required in the following. Now the inverse 4 X 1 is simply the diagonal matrix with 
entries 1/Xi(u) and, given eq. (6.8), Xa b is also diagonal with 

X A B = diag (x^^^u),^),^),^),^))^ 

f 1 1 1 1 1 1 \ 

Now we review the general structure of the Green's functions in the boundary 
theory, from the discussion in [7] . Together current conservation and spatial rotational 
invariance - Lorentz invariance is lost with T ^ - dictate the form of the retarded 
Green's functions as 

= V 7 ? (Pi K T (co, q) + I*, K L (u, q)) . (A.4) 

where we use the notation: q M = (u,q x ,q y ), q 2 = \{q x ) 2 + (q v ) 2 } 1 ^ 2 and q 2 = q 2 — tu 2 . 
Further, P^ u and P^ v are orthogonal projection operators defined by 

Pit = = P ti = P it , P t j = dij — , P^ u = I rj^ — T^p- \ — P^ u , (A.5) 

with i, j denoting spatial indices while \i , v run over both space and time. If, for 
simplicity, we choose q M = (u,q,0), then we have 

Q yy (u,q) = Vq 2 -u 2 K T (u,q), g tt (u,q) = - J? K L (u,q). (A.6) 

\/ q A — u z 

Of course, this general structure applies for both boundary theories, that is, both for 
the theory dual to the vector potential A a and that dual to B a . Our notation will be 
that the above expressions refer to the theory dual to A a while Q^, K T and K L are 
the corresponding expressions for the boundary currents dual to B a . 

The first step in the holographic calculation of the Green's functions is to solve 
the bulk equations of motion. Hence we begin as in section 4 by taking a plane-wave 
ansatz (4.1) for A a and B a . Further, we choose q M = (oj, q, 0) and work in radial gauge 
with A u (u, q) = = B u (u, q). With these choices and the background metric (2.4), 

4 We assume that the functions Xi remain finite and positive throughout u £ [0, 1] in order that 
Xa b is invertible and the bulk propagators for the gauge potential are well-behaved there. 
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the A a equations of motion become: 



A' t + q ^A' x = (A.7) 



u X 



A'; + §A' t -^ q -^ (gA t + uA x ) = ( A.8) 

^+[^+ f -?iK + K^(QA t + ujA x ) = (A.9) 



5 



X 5 1 fj x r 2 f 2 X { 
X> , /'\ „, L*u 2 X 2 -q 2 fX 4 



A « + {t + 7) A »-^^px~ A » = ° (A ' 10) 

where we recall that / = 1 — w 3 . For the EM dual gauge theory, the equations of motion 
are given by simply replacing A a — > B a and Xi — > Xi in the expressions above. 

In general, there are two independent physical modes for the four-dimensional 
bulk gauge field. Above, we see that A y decouples in eq. (A. 10) to provide one of 
these modes, while A t and A x are coupled in the remaining equations. Of course, the 
analogous results apply to B a in the EM dual theory. Now explicitly writing out the 
duality relations (6.15) in the present case, we find 



fa 



9± r _^L r p 9 * r of 

X 1 L 2 yu ' ty ' X 2 L 2 

A^r tp - -AL 

X 3 r xy ' xy ~ X 4 L 

2 j 2 2 j 2 

9a ^ n Tp 9 a ^ 



G xy , F xy — —^-—Gtu, (A. 11) 



X 5 r f X 6 r f 

Hence, at a schematic level, EM duality exchanges the A y mode for that in B t , x and 
similarly the A ttX and B y are exchanged. Given the holographic relationship between 
the bulk and boundary theories, we expect that there are connections between the 
Green's functions, Q^ v and Q^, generalizing those found in [7]. However, given the 
previous observation, more specifically, Q yy should be related to Q t t (as well as Q xx and 
Qtx) and similarly Q tu to Q yy . 

To develop these connections in detail, we must extend the holographic calculation 
of the Green's functions given in section 4 to include the mixing between A t and A x , 
noted above. First, we solve the equations of motion (A.7)-(A.10) for A^ with infalling 
boundary conditions at the horizon and asymptotic boundary conditions: lim A^ = A®. 
To account for mixing between different components of the gauge potential, we may 
write [7]: A^iu) = M M "(tt)A°. Now, substituting the solutions into the action (3.1) 
and integrating by parts leaves an surface term at the asymptotic boundary, which 
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generalizes that given in eq. (4.12) 

2ttT 



h = 2 ~mS d " x [ x 3M-^4a:-i 6 v;]^ . (a. 12) 

After Fourier transforming in the boundary directions, we extract the desired Green's 
functions as 



Gxx{u,q) 

Gyy(u,q) 



4vrT 



4ttT 
2tvT ~ 



X 3 (0) 



SA' t (u) 



SA9 



«->0 



X 5 (0) 



SA'Ju) 



6A1 



3gl 



AttT 



XJO) 



5AI 
SA'Ju) 



5A° 



t J u->0 



6A° 



(A.13) 
(A. 14) 
(A.15) 
(A.16) 



Here we have used that the equations of motion (A.7-A.10) only mix A t and A x . One 
may also easily verify that eq. (A.16) reduces to the expression in eq. (4.16) when 
Xq(0) = 1, as in the main text. 

Next consider the Green's functions Q m . Assume that we have A y (u) = ip(u)Ay 
where ip(u) is a solution of eq. (A. 10) satisfying the appropriate boundary conditions. 
In particularly, the asymptotic normalization is ip(u = 0) = 1. Then from (A.16), we 
have 

A-kT 

g w (u,q) = -- T M0)^(P)- (A- 17) 

Given that EM duality exchanges A y with B t)X , we now look for a relation between 
this result and that for Q tt - From the expression for F xy in eq. (A. 11), we find B' t {u) oc 
X±(u) A y {u) and so X 4 (w) ip(u) provides a solution of the equations of motion for B' t {u) 
in the EM dual theory. While it is clear that the required infalling boundary condition 
is satisfied at the horizon with ip(u), we must expect that the normalization has to 
be adjusted in order to satisfy the desired asymptotic boundary condition. Hence 
we introduce a new constant C\ setting B[{u) = C\ Xju) ip(u). In order to fix this 
constant, we consider the analog of eq. (A. 8) in the EM dual theory and take the limit 
u — > to find 



Ci 



LU{qB° + uBl 



where deriving this expression uses X\ 
counterpart of eq. (A.13) yields 



2 X 6 (0)^(0) 

I/Xq and X 3 



(A.18) 

1/X 4 . Now the EM dual 



- 34 - 



p 4ttT - . , SB'Ju) 

Gtt(oJ,q) = ^ Y X 3 (0) 



6B? 



M->0 



3^ 2 g 2 
4ttT X 6 (0)^'(0) ' 



(A.19) 



Here we have used the relations: = 1/(74 and tq/L 2 = 4ttT/3. Hence, combining 
eqs. (A. 17) and (A.19), we find 



Gyy{u,q)Gtt{u,q) = ~q 2 ■ 
Further, using eq. (A. 6), this relation can be written as 

K T (u,q)K L (u,q) = l. 



(A.20) 



(A.21) 



Now it is clear that the EM dual version of the above discussion would follow 
through without change. That is, we would begin by constructing an expression for 
G yy analogous to eq. (A. 17) and then the counterpart of eq. (A.19) for Gtt- The final 
result emerging from these results would then be 



K T (u,q)K\u,q) = l. 



(A.22) 



To close our discussion, we comment on more general cases where X contains off- 
diagonal terms. To begin, let us write the most general tensor which is consistent with 
rotational symmetry in the xy-plane: 



X A B 



u 



L 2 7 

¥ \ 



XAu) 



— f z 2 {u) 



X 3 (u) 



ra 

X 



rof 
L 2 



«) 

U) 



rof 



zi(u) 



X, 



u 



(A.23) 



where we are using the notation introduced in eq. (A.l), as well as the background 
metric (2.4). Note the pre-factors in the off-diagonal terms reflect the tensor structure 
of X a b cd , which is slightly obscure in this notation, e.g., X5 1 = g xx g uu g u g xx X\ 5 = 
— L 4 /(ro/) 2 A4 5 . Now, as noted above, rotational invariance imposes two relations 
on the diagonal entries, i.e., X 2 = X\ and X$ = X 5 . However, as shown above, 
this symmetry is remarkably restrictive on the off-diagonal components as well and 
our general tensor (A.23) only contains three independent terms amongst all of the 
possible entries. Now, if we further demand that this background tensor preserves 
parity, we must in fact set z 2 (u) = = z 3 (u) and we are left with only one function 
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Z\ (it) determining all of the allowed off-diagonal components. Note that these remaining 
off-diagonal terms preserve parity but violate time-reversal invariance. 5 

If we restrict ourselves to the parity invariant case, it is straightforward to generalize 
our previous discussion to accommodate the general X above (with Z2 = = 23). 
Although the intermediate expressions are somewhat more involved, we find that the 
final Green's functions still satisfy eqs. (A. 21) and (A. 22). 

Note that parity invariance was implicit in the decomposition of the Green's func- 
tions in eq. (A. 4). If parity violating terms were allowed there would be an additional 
contribution of the form. 

&g^ = ie^ a q a K P (oo,q). (A.24) 

Hence the present analysis must be revised to accommodate these parity violating 
terms. Our expectation is that particle-vortex duality still provides relations between 
the three functions K T , K L and K p , describing the Green's functions of the two dual 
theories. A preliminary examination of the equations of motion and the EM duality 
relations suggests that, in this general case, K T , K L , K p and their dual counterparts 
should satisfy three relations. However, the details of this interesting case are left as 
an open problem for future work. 
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